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Abstract Being the limits of copulas of componentwise maxima in independent 
random samples, extreme- value copulas can be considered to provide appropriate 
models for the dependence structure between rare events. Extreme-value copulas 
not only arise naturally in the domain of extreme-value theory, they can also be a 
convenient choice to model general positive dependence structures. The aim of this 
survey is to present the reader with the state-of-the-art in dependence modeling via 
extreme-value copulas. Both probabilistic and statistical issues are reviewed, in a 
nonparametric as well as a parametric context. 
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1 Introduction 

In various domains, as for example finance, insurance or environmental science, 
joint extreme events can have a serious impact and therefore need careful modeling. 
Think for instance of daily water levels at two different locations in a lake during a 
year. Calculation of the probability that there is a flood exceeding a certain bench- 
mark requires knowledge of the joint distribution of maximal heights during the 
forecasting period. This is a typical field of application for extreme-value theory. 
In such situations, extreme-value copulas can be considered to provide appropriate 
models for the dependence structure between exceptional events. 
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One of the first applications of bivariate extreme-value analysis must be due to 
Gumbel and Goldstein 14711 . who analyze the maximal annual discharges of the Oc- 
mulgee River in Georgia at two different stations, a dataset that has been taken up 
again in M56I1 . The joint behavior of extreme returns in the foreign exchange rate 
market is investigated in J87ll . whereas the comovement of equity markets charac- 
terized by high volatility levels is studied in l65tl . An application in the insurance 
domain can be found in Jgt] . 

Extreme-value copulas not only arise naturally in the domain of extreme events, 
but they can also be a convenient choice to model data with positive dependence. An 
advantage with respect to the much more popular class of Archimedean copulas, for 
instance, is that they are not symmetric. Incidentally, a hybrid class containing both 
the Archimedean and the extreme-value copulas as a special case are the Archimax 
copulas |@]. 

The aim of this survey is to present the reader with the state-of-the-art in depen- 
dence modeling via extreme-value copulas. Definition, origin, and basic properties 
of extreme-value copulas are presented in Section [2] A number of useful and pop- 
ular parametric families are reviewed in SectionfJ] Section @]provides a discussion 
of the most important dependence coefficients associated to extreme-value copulas. 
An overview of parametric and nonparametric inference methods for extreme-value 
copulas is given in Section[5] Finally, some further topics and pointers to the litera- 
ture are gathered in Section[6] 



2 Foundations 

Let Xt = (Xfl , . . . jXtf), i G { 1 , . . . , «}, be a sample of independent and identically 
distributed (iid) random vectors with common distribution function F, margins 
F\,... ,Fa, and copula Cf- For convenience, assume F is continuous. Consider the 
vector of componentwise maxima: 

n 

M n = (M„ il; ...,M M ), where M nJ = \J X ih (1) 

/=i 

with 'V' denoting maximum. Since the joint and marginal distribution functions of 
M„ are given by F" and F",..., F? respectively, it follows that the copula, C„, of M„ 
is given by 

C„(u u ...,u d ) =C F (u l /'\...,u l J")" 7 (u u ...,u d ) G [0,l} d - 

The family of extreme-value copulas arises as the limits of these copulas C„ as the 
sample size n tends to infinity. 

Definition 2.1. A copula C is called an extreme-value copula if there exists a copula 
Cf such that 
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CF(u l /",...,uJ n ) n -^C{u\,...,u d ) (n-»-oo) (2) 

for all (u\ , . . . , iifi) £ [0, l] d . The copula Cf is said to be in the domain of attraction 
ofC. 



Historically, this construction dates back at least to 11161 13011 . The representation 
of extreme-value copulas can be simplified using the concept of max-stability. 

Definition 2.2. A (i-variate copula C is max-stable if it satisfies the relationship 

C{u u ...,u d ) =C{u{ ,...,uj ) (3) 

for every integer m 1 and all {u\ , . . . , u d ) £ [0, l] d . 

From the previous definitions, it is trivial to see that a max-stable copula is in 
its own domain of attraction and thus must be itself an extreme-value copula. The 
converse is true as well. 

Theorem 2.1. A copula is an extreme-value copula if and only if it is max-stable. 



The proof of Theorem 12. II is standard: for fixed integer m ^ 1 and for n = mk, 
write 

C„{ui,...,u d ) =C k {u{ ,...,uj ) . 

Let k tend to infinity on both sides of the previous display to get OJ. 

By definition, the family of extreme-value copulas coincides with the set of cop- 
ulas of extreme-value distributions, that is, the class of limit distributions with non- 
degenerate margins of 

M„ } i-b„ t i M„ :d - b n4 

3 • • • 5 

with M n j as in (Q]i, centering constants b n j £ K and scaling constants a n j > 0. Rep- 
resentations of extreme-value distributions then yield representations of extreme- 
value copulas. Let A d ^\ — {{w\ , . . -,w d ) £ {0,°°) d : L; W; = 1} be the unit simplex 
in W 1 ; see FigureQ] The following theorem is adapted from $T]\\ . which is based in 
turn on Jill . 

Theorem 2.2. A d-variate copula C is an extreme-value copula if and only if there 
exists a finite Borel measure H on A d -h called spectral measure, such that 

C(m,...,u d ) =exp(—£(-logwi,...,-logM d )), (u u ...,u d ) £ (0,1] (/ , 

where the tail dependence function £ : [0,°°) d — > [0,°°) is given by 
d 

£(x u ...,x d )= \/(wjXj)dH(w h ...,w d ), (xi,...,x d )e[0,°o) . ( 4 ) 
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e 3 = (0,0,1) 
A 2 

62 = (0,1,0) 

eipji',0,0) 




c 2 = (0, 1) 

^_ei = (l,0) 







Fig. 1 On the left side is represented in R 3 , which is equivalent to the representation in K 2 on 
the right side. 



The spectral measure H is arbitrary except for the d moment constraints 

/ w j dH(w 1 ,...,w d ) = l, je{l,...,d}. (5) 

The d moment constraints on H in stem from the requirement that the margins 
of C be standard uniform. They imply that H{A d _\) = d. 

By a linear expansion of the logarithm and the exponential function, the domain- 
of-attraction equation (f2]i is equivalent to 

]har 1 (l-C F (l-tx 1 ,...,l-tx d )) = -logC(e~ Xl ,...,e~ Xi ) 

fj.0 

= £(x l ,...,x d ) (6) 

for all (jci, . . . ,x d ) e [0,00)^; see for instance fl9ll . The tail dependence function I in 
(0]i is convex, homogeneous of order one, that is l(cx\ cx d ) = ci{x\ ,x d ) for 

c> 0, and satisfies max(xi ,...,x d ) ^ l(x\ ,x d ) ^x\A \-x d for all (x\ , . . . ,x d ) 6 

[0,°°) d . By homogeneity, it is characterized by the Pickands dependence function 
A : A ( [-i — > [1/d, 1], which is simply the restriction of t to the unit simplex: 

i(x 1 ,...,x d ) = (*H \-x d )A(wi,...,w d ) where w,- = , 

x\ H hXd 

for (xi, . . . ,x^) € [0,oo) </ \ {0}. The extreme-value copula C can be expressed in 
terms of A via 

C(m,...,u d )= ex.pl I £ log itj ] A ( rf l0g , ■ ■ ■ , V , ; l0g - Mt/ ) [• 
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The function A is convex as well and satisfies max(wi , . . . , w^) ^ A(wi , ...,Wd) ^ 1 
for all (w\ ,Wd) £ A^-i- However, these properties do not characterize the class 
of Pickands dependence functions unless d = 2, see for instance the counterexample 
on p. 257 in Jj]. 

In the bivariate case, we identify the unit simplex A\ = {(1 — t,t) : t £ [0, 1]} in 
R 2 with the interval [0, 1]. 

Theorem 2.3. A bivariate copula C is an extreme-value copula if and only if 

C( M ,v) = ( M v) A(log(v)/log( " v)) , («,v)G(0,l] 2 \{(l,l)}, (7) 

where A : [0, 1] — ► [1 /2, 1] is convex and satisfies t V (1 — /) ^ A(f) ^ 1 for all t £ 
[0,1]. 

It is worth stressing that in the bivariate case, any function A satisfying the two 
constraints from Theorem l2 . 3 1 corresponds to an extreme-value copula. These func- 
tions lie in the shaded area of Figure |2j in particular, A(0) = A ( 1 ) = 1. 



A(t) 



Fig. 2 A typical Pickands 
dependence function A 
together with the region 
fV(l-f) <A(f) =S 1 in 
Theorem l2.3l 



The upper and lower bounds for A have special meanings: the upper bound 
A(f) = 1 corresponds to independence, C(u,v) = uv, whereas the lower bound 
A(f) = t V (1 — t) corresponds to perfect dependence (comonotonicity) C(u,v) = 
u A v. In general, the inequality A (f) ^ 1 implies C(u,v) uv, that is, extreme-value 
copulas are necessarily positive quadrant dependent. 



3 Parametric models 

By Theorems |2.2| and [2.3l the class of extreme-value copulas is infinite-dimensional. 
Parametric submodels can be constructed in a number of ways: by calculating the 
limit £ in © for a given initial copula CV; by specifying a spectral measure H; in 
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dimension d = 2, by constructing a Pickands dependence function A. In this section, 
we employ the first of these methods to introduce some of the more popular families. 
For more extensive overviews, see e.g. 15 , 62 1 . 



3.1 Logistic model or Gumbel-Hougaard copula 



Consider the Archimedean copula 

C$(ui,...,u d ) = (j>^((j>(ui) + --- + (j>(u d )), (ui,...,u d ) £ [0,l] d (8) 

with generator (j> : [0, 1] — > [0,°°] and inverse (j>^(t) = inf{« G [0, 1] : <j)(u) ^ f}; the 
function should be strictly decreasing and convex and satisfy 0(1) =0, and 0^ 
should be cZ-monotone on (0,°°), see 116911 . 
If the following limit exists, 



40 0(1-5) 1 



(9) 



then the domain-of-attraction condition © is verified for Cf equal to Cq, the tail 
dependence function being 



Ui VL-.-.A, if0=oc, 



(10) 



for (jci, . . . ,x f /) G [0,oo) d ; see Jl, 11 1. The range [1,°°] for the parameter 9 in © is 
not an assumption but rather a consequence of the properties of (p. The parameter 6 
measures the degree of dependence, ranging from independence (9 = 1) to complete 
dependence (9 = °°). 

The extreme-value copula associated to £ in ( fTOl is 

C( Ml ,...,M rf )=exp{-((-log Ml ) e + -.. + (-logM (/ ) e ) 1/e }, 



known as the Gumbel-Hougaard or logistic copula. Dating back to Gumbel J44 , 
4511 . it is (one of) the oldest multivariate extreme-value models. It was discovered 
independently in survival analysis fl4jl52Tl . It hap pen s to be the only copula that is 
at the same time Archimedean and extreme-value [35]. 



The bivariate asymmetric logistic model introduced in [91 ] adds further flexibil- 
ity to the basic logistic model. Multivariate extensions of the asymmetric logistic 
model were studied already in 16811 and later in 11131 15811 . These distributions can 
be generated via mixtures of certain extreme-value distributions over stable distri- 
butions, a representation that yields large possibilities for modelling that have yet 
begun to be explored 1281 19211 . 
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3.2 Negative logistic model or Galambos copula 

Let C(j, be the survival copula of the Archimedean copula Ca in ([8]). Specifically, 
if Ca is the distribution function of the random vector (U\ , t/ f /), then Ca is the 
distribution function of the random vector (1 — Ui, . . . , 1 —U c j). If the following limit 
exists, 



e 



lim- 



G[0, 



(11) 



slO S(j)'(s) 

then the domain-of-attraction condition (O is verified for CV equal to Ca, the tail 
dependence function being 



e(x h ...,x d ) = < 



X] 



Xd- 



E (-i) |7, (Ete/*, 



-e\-i/fl 



'c{i rf} 

1/1^2 



xi V • • • Vij 



if e = o, 

if < <°c, 
if 9 =oo, 



for (xi, . . . ,Xrf) G [0,o°) f/ ; see SEIll- In case < < oo, the sum is over all subsets 
/of {l,...,d} of cardinality |/| at least 2. The amount of dependence ranges from 
independence (9 = 0) to complete dependence (9 = °°). 

The resulting extreme-value copula is known as the Galambos or negative lo- 
gistic copula, dating back to I29I1 . Asymmetric extensions have been proposed in 

EH. 



3.3 Hiisler-Reiss model 

For the bivariate normal distribution with correlation coefficient p smaller than one, 



it is known since II 8411 that the marginal maxima M„,\ and M, h 2 are asymptotically 
independent, that is, the domain-of-attraction condition (O holds with limit cop- 
ula C(u,v) = uv. However, for p close to one, better approximations to the copula 
of M n ,\ and M„ 9 arise within a somewhat different asymptotic framework. More 
precisely, as in |543, consider the situation where the correlation coefficient p as- 
sociated to the bivariate Gaussian copula C p is allowed to change with the sample 
size, p = p„, in such a way that p„ — > 1 as n — > 00. If 

(l-p„)log«^A 2 e [0,+-] (n -><»), 

then one can show that 

C A ,(« 1 /» v i/»)»_>c A («,v) (»-^~). («,v)e[0,l] 2 , 



where the Hiisler-Reiss copula Ca is the bivariate extreme- value copula with Pickands 
dependence function 
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A{w) = (l-w)<P A + — log +w<t> U + — log- 

\ 2/i w J \ 2A 1 —w 

for w £ [0, 1], with <£> representing the standard normal cumulative distribution func- 
tion. The parameter A measures the degree of dependence, going from independence 
(A = °°) to complete dependence (A = 0). 



3.4 The t-EV copula 

In financial applications, the f-copula is sometimes preferred over the Gaussian cop- 
ula because of the larger weight it assigns to the tails. The bivariate f-copula with 
V > degrees of freedom and correlation parameter p £ (— 1 , 1 ) is the copula of the 
bivariate f-distribution with the same parameters and is given by 



i r(Jf + i)/ x 'p- l x\-^ 2+l J 



where t v represents the distribution function of the univariate f-distribution with v 
degrees of freedom and P represents the 2x2 correlation matrix with off-diagonal 



element p. In III 811 . it is shown that C v .p is in the domain of attraction of the bivariate 



extreme-value copula Ca with Pickands dependence function 

A(w) =wf v +i(z w ) + (l -w)f v +i(zi- w ), 
wherez vl . = (l + v) 1 / 2 [{w/(l-w)} 1 / v -p](l-p 2 )- 1 / 2 , we [0,1]. (12) 



This extreme-value copula was coined the t-EV copula. Building upon results in |3S_ 
511. exactly the same extreme-value attractor is found in J2l for the more general 



class of (meta-)elliptical distributions whose generator has a regularly varying tail. 



4 Dependence coefficients 

Let (U,V) be a bivariate random vector with distribution function C, a bivariate 
extreme-value copula with Pickands dependence function A as in @. As mentioned 
already, the inequality A ^ 1 implies that C(u,v) ^ uv for all (w, v) G [0, l] 2 , that is, 
C is positive quadrant dependent. In fact, in 114211 it was shown that extreme-value 
copulas are monotone regression dependent, that is, the conditional distribution of U 
given V — v is stochastically increasing in v and vice versa; see also Theorem 5.2.10 



in ||78ll. 

In particular, all measures of dependence of C such as Kendall's T or Spearman's 
Ps must be nonnegative. The latter two can be expressed in terms of A via 
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t = 4 /|, iP c <"'" )dc( "^ ) - 1 = r^ dA '(')' 

ft = 12 /i,p"' ,dC( "'' ,) - 3 = 12 i'(rpiMP d '- 3 - 

The Stieltjes integrator dA'(f) is well-defined since A is a convex function on [0, 1]; 
if the dependence function A is twice differentiable, it can be replaced by A"(t) dt . 
For a proof of the identities above, see for instance ll53ll . where it is shown that T 
and ps satisfy — 1 + ^/T+3r p s < min (|t,2t — t 2 ), a pair of inequalities first 
conjectured in 15511 . 

The Kendall distribution function associated to a general bivariate copula C is 
defined as the distribution function of the random variable C(U,V), that is, 

K(w)=Y[C(U,V)^w\, we [0,1]. 

The reference to Kendall stems from the link with Kendall's T, which is given by T = 
4 E[C(U, V)] — 1. For bivariate Archimedean copulas, for instance, the function K 
not only identifies the copula [36], convergence of Archimedean copulas is actually 
equivalent to weak convergence of their Kendall distribution functions ifioll . For 
bivariate extreme-value copulas, the function K takes the remarkably simple form 

K(w)=w-(l-T)w\ogw, we [0,1], (13) 

as shown in fiiil . In fact, in that paper the conjecture was formulated that if the 
Kendall distribution function of a bivariate copula is given by 1131 . then C is a 
bivariate extreme-value copula, a conjecture which to the best of our knowledge 
still stands. In the same paper, equation ( fT3l was used to formulate a test that a 



copula belongs to the family of extreme-value copulas; see also [38]. 

In the context of extremes, it is natural to study the coefficient of upper tail de- 
pendence. For a bivariate copula Cp in the domain of attraction of an extreme-value 
copula with tail dependence function £ and Pickands dependence function A, we 
find 

X v = limP(t/ > u I V > u) = \\mr x ( 2t- 1 +C(1 - 1, 1 -t)) 
«Ti fj.0 

= 2-£(l,l)=2(l-A(l/2))e[0,l]. 

Graphically this quantity can be represented as the length between the upper bound- 
ary and the curve of the Pickands dependence function evaluated in the mid-point 
1/2, see Figure [3] The coefficient Ay ranges from (A = 1, independence) to 1 



(complete dependence). Multivariate extensions are proposed in [64]. 

The related quantity ^(1,1) = 2A(1 /2) is called the extremal coefficient in II86I1 
For a bivariate extreme-value copula, we find 

P(U^u,V s^m) = m 2a(1/2) , «e [0,1], 
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Fig. 3 The coefficient of I 



upper tail dependence Xy is 
equal to twice the length of 
the double arrow in the upper 
part of the graph. 




t 



so that 2A(1 /2) g [1,2] can be thought of as the (fractional) number of i ndep endent 
components in the copula. Multivariate extensions have been studied in 18011 . 
For the lower tail dependence coefficient, the situation is trivial: 

X L = BmP(CT <u\V< u) = hW^/2)-D = J° if A(1/2) > 1/2 ' 
»|0 40 (1 ifA(l/2) = 1/2. 

In words, except for the case of perfect dependence, A(l/2) = 1/2, extreme-value 
copulas have asymptotically independent lower tails. 



5 Estimation 

Let X, = (Xj\ , . . . ,Xj t i), i <E {1, . . . ,«}, be a random sample from a (continuous) dis- 
tribution F with margins F\ , . . . , and extreme-value copula C: 

F(xi, . . . ,x d ) = C(Fi (xi ), . . .,F d (xd)) , 

and C as in Theorem |2.2| The problem considered here is statistical inference on C, 
or equivalently, on its Pickands dependence function A. A number of situations may 
arise, according to whether the extreme-value copula C is completely unknown or is 
assumed to belong to a parametric family. In addition, the margins may be supposed 
to be known, parametrically modelled, or completely unknown. 
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5.1 Parametric estimation 



Assume that the extreme-value copula C belongs to a parametric family (Cg : 9 € 0) 
with © C M. p ; for instance, one of the families described in Section[3] Inference on 
C then reduces to inference on the parameter vector 9. The usual way to proceed is 
by maximum likelihood. The likelihood is to be constructed from the copula density 

d d 

c e (ui,...,u d ) = 3 5 — C g (ui,...,u d ), (ui,...,u d ) G (0,l) d . 

au\ ■ ■ ■ ou d 

In order for this density to exist and to be continuous, the spectral measure H should 
be absolute continuous with continuous Radon-Nikodym derivative on all 2^—1 
faces of the unit simplex with respect to the Hausdorff measure of the appropriate 



dimension H13H . In dimension d = 2, the Pickands dependence function A : [0, 1] — > 
[1/2,1] should be twice continuously differentiable on (0,1), or equivalently, the 
spectral measure H should have a continuous density on (0, 1) (after identification 
of the unit simplex in M 2 with the unit interval). 

In case the margins are unknown, they may be estimated by the (properly 
rescaled) empirical distribution functions 

F nj {x) = -^-f i I(X ij ^x), xeR, j €{!,.. .,d}. (14) 
n + 1 (=1 

(The denominator is n + 1 rather than n in order to avoid boundary effects in the 
pseudo-loglikelihood below.) Estimation of 9 then proceeds by maximizing the 
pseudo-loglikelihood 

n 

£ logc e (F nX (X a ), . . .,F nd (X id )) , 



see 13411 . The resulting estimator is consistent and asymptotically normal, and its 
asymptotic variance can be estimated consistently. 

If the margins are modelled parametrically as well, a fully parametric model for 
the joint distribution F arises, and the parameter vector of F may be estimated by 
ordinary maximum likelihood. An explicit expression for the 5x5 Fisher informa- 
tion matrix for the bivariate distribution with Weibull margins and Gumbel copula is 



calculated in 117 1H . A multivariate extension and with arbitrary generalized extreme 



value margins is presented in 18311 . 

Special attention to the boundary case of independence is given in 19111 . In this 
case, the dependence parameter lies on the boundary of the parameter set and the 
Fisher information matrix is singular, implying the normal assumptions for validity 
of the likelihood method are no longer valid. 

A robustified version of the maximum likelihood estimator is introduced in (2 



The effects of misspecification of the dependence structure are studied in [21]. 
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5.2 Nonparametric estimation 



For simplicity, we restrict attention here to the bivariate case. For multivariate ex- 
tensions, see BUB. 

Let (X\ , Y\ ) , . . . , (X n , Y„ ) be an independent random sample from a bivariate dis- 
tribution F with extreme-value copula C and Pickands dependence function A. As- 
sume for the moment that the marginal distribution functions F\ and Fj are known 
and put Uj = Fi (X t ) and V, = F 2 (Jt ) and 5, = - log t/,- and 7} = - log V, . Note that S, 
and 7} are standard exponential random variables. For t 6 [0, 1], put 



Si T, 



1-f' t 



with the obvious conventions for division by zero. A characterizing property of 
extreme-value copulas is that the distribution of 4/(0 is exponential as well, now 
with mean 1/A(f): for x > 0, 

P[&(f) > x] = P[U, < e-^Wi < e-' x ] 

= C(e-( 1 - t >,e- tx )=e- xA ®. (15) 



This fact leads straightforwardly to the original Pickands estimator 17711 : 

A major drawback of this estimator is that it does not verify any of the constraints 
imposed on the family of the Pickands dependence functions in Theorem [23] 

Besides establishing the asymptotic properties of the original Pickands estimator, 
Deheuvels lfl7ll proposed an improvement of the Pickands estimator that at least 
verifies the endpoint constraints A (0) = A(l) = 1: 

3 4 = it6(')-'iE6(i)-a-«)E6(0) + 1 . (17) 



As shown in 18211 . the weights (1 — t) and t in de Deheuvels estimator (QjJ) can be 
understood as pragmatic choices that could be replaced by suitable weight functions 
j3 1 (f)andj3 2 (f): 

- £ $(/) - p^o l - £ £(i) - ftw t &(°) + L < 18 > 

11 "™ 11 "™ "™ 



The linearity of the right-hand side of ( fT7] i in 4/(0 suggests to estimate the variance- 
minimizing weight functions via a linear regression of £,•(/) upon 4/(0) and 4/(1): 



m = po( t ) +ft(o {4,(o) - i}+&(o {4,(1) - 1} +m- 
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The estimated intercept j3o(f) corresponds to the minimum- variance estimator for 
1 /A(t ) in the class of estimators (fT~8T >. 

In the same spirit, Hall and Tajvidi l50ll proposed another approach to improve 
the small-sample properties of the Pickands estimator at the boundary points. For 
all / G [0, 1] and i £ {1, . . .,»}, define 



|,-(f) =min( 



1-r' t 



with 



i(Si+-- + S n y l( Tl + ... + T n ) 

The estimator presented in f50ll is given by 

1 1 " - 

ThtTT = ~ / , v(f) ■ 
A (t) »g 

Not only does the estimator's construction guarantee that the endpoint conditions are 
verified, in addition it always verifies the constraint A HT (t) 1 V (1 — t). Among 
the three nonparametric estimators mentioned so far, the Hall-Tajvidi estimator typ- 
ically has the smallest asymptotic variance. 

A different starting point was chosen by Caperaa, Fougeres and Genest nfl: they 
showed that the distribution function of the random variable Z, = log(I//)/log(I/,-V/) 
is given by 

P(ZKz)=z + z(l-z)^, (Kz<l, 
A(z) 

where A' denotes the right-hand derivative of A. Solving the resulting differential 
equation for A and replacing unknown quantities by their sample versions yields the 
CFG-estimator. In 18211 however, it was shown that the estimator admits the simpler 
representation 

logA CFG (f) = -±J>g&(f) - (1 -t)f>g&(0)-f f>g£(l) (19) 

n i=\ i=i i=i 

for t € [0, 1]. This expression can be seen as a sample version of 

E[-log6(0]=logA(*) + r, re [0,1], 

a relation which follows from (fl3T i: note that the Euler-Mascheroni constant y = 
0.5772 ... is equal to the mean of the standard Gumbel distribution. Again, the 
weights (1 — f) and t in (fT9l can be replaced by variance-minimizing weight func- 



tions that are to be estimated from the data 114 1L 18211 . The CFG-estimator is con 



by v 
EI 



sistent and asymptotically normal as well, and simulations indicate that it typically 
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performs better than the Pickands estimator and the variants by Deheuvels and Hall- 
Tajvidi. 

Theoretical results for extreme values in the case of unknown margins are quite 



recent. To some extent, Jimenez, Villa-Deharce and Flores H56H were the first to 
present an in-depth treatment of this situation. However, their main theorem on uni- 
form consistency is established under conditions that are unnecessarily restrictive. 
In ll37tl . asymptotic results were established under much weaker conditions. The es- 
timators are the same as the ones presented above, the only difference being that 
Uj = F\ (X{) and V; = FjiYf) are replaced by 

with F n j as in <fT~4t . Observe that the resulting estimators are entirely rank-based. 
Contrary to the case of known margins, the endpoint-corrections are irrelevant in 
the sense that they do not show up in the asymptotic distribution. Again, the CFG- 
estimator has the smallest asymptotic variance most of the time. 

The previous estimators do typically not fulfill the shape constraints on A as 
given in Theorem l2.3l A natural way to enforce these constraints is by mo difying a 



lity i 

pilot estimate A into the convex minorant of (A(t) V (1 — t) Vt) A 1, see Ill7l.l56l.l77 l. 
It can be shown that this transformation cannot cause the L°° error of the estimator 
to increase. A different way to impose the shape constraints is by constrained spline 
smoothing [01 or by constrained kernel estimation of the derivative of A lHol . 

The L 2 -viewpoint was chosen in l26ll . The set s/ of Pickands dependence func- 
tions being a closed and convex subset of the space L 2 ([0, 1], dx), it is possible to 



find for a pilot estimate A a Pickands dependence function Aea' that minimizes 
the L 2 -distance Jq(A—A) 2 . By general properties of orthogonal projections, the 
L 2 -error of the projected estimator cannot increase. 

Finally, a nonparametric Bayesian approach has been proposed by Guillotte and 
Perron 114311 . Driven by a nonparametric likelihood, their methodology yields an es- 
timator with good properties: its estimation error is typically small, it automatically 
verifies the shape constraints, and it blends naturally with parametric likelihood 
methods for the margins. 



6 Further reading 



About the first monograph to treat multivariate extreme-value dependence is the one 
by Galambos J30ll . with a major update in the second edition |31|. Extreme-value 
copulas are treated extensively in the monographs J5 , 62 ] and briefly in the 2006 edi- 
tion of Nelsen's book lf7(ill . The regular- variation approach to multivariate extremes 
is emphasized in the books by Resnick 178i 17911 and de Haan and Ferreira 114 8H . A 
highly readable introduction to extreme-value analysis is the book by Coles 11211 . 
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The first rep resentations of bivariate extreme-value distributions are due to 
Finkelstein J27D, Tiago de Oliveira H, Geffroy JHO and Sibuya JH. Inciden- 



tally, the 1959 paper by Geffroy appeared in the same issue as the famous paper by 
Sklar 1 85]. The equivalence of all these representations was shown in Gumbel j46ll ; 
see also the more recent paper by Obretenov [72]. However, their representations of 
multivariate extreme value distributions have not enjoyed the same success as the 
one proposed by Pickands 17711 . The domain of attraction condition seems to have 
been formulated for the first time by Berman 0], his standardization being to the 
standard exponential distribution rather than the uniform one. 

A particular class of extreme-value copulas arises if the spectral measure H in 
Theorem l2.2l is discrete. In that case, the stable tail dependence function I and the 
Pickands dependence function A are piecewise linear. In general, such distributions 
arise from max-linear combinations of independent random variables llioll . An early 
example of such a distribution is the bivariate model studied by Tiago de Oliveira in 
ll74l[75 . 761. which has a spectral measure with exactly two atoms; see also 1 24 ] . The 
bivariate distribution of Marshall and Olkin |67ll corresponds to a spectral measure 
with exactly three atoms, jO, 1 /2, 1} (after identification of the unit simplex in R 2 
with the unit interval); see 16611 for a multivariate extension. 

Even more challenging than the estimation problem considered in Section [5] is 
when the random sample comes from a distribution which is merely in the domain 
of attraction of a multivariate extreme-value distribution. See for instance |5 j I 1 3l |49t 
EUS EH El for some (semi-)parametric approaches and |fl 151 Izzj \2S,m\\ 
for some nonparametric ones. 

For an overview of software related to extreme value analysis, see JUt]. Par- 
ticularly useful are the R packages evd HU, which provides algorithms for the 
computation, simulation |9()| and estimation of certain univariate and multivariate 
extreme-value distributions, as well as the more general copula package 193 1. 
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